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We propose an extension of natural inflation, where the inflaton potential is a general periodic 
function. Specifically, we study elliptic inflation where the inflaton potential is given by Jacobi 
elliptic functions, Jacobi theta functions or the Dedekind eta function, which appear in gauge and 
Yukawa couplings in the string theories compactifled on toroidal backgrounds. We show that in the 
first two cases the predicted values of the spectral index and the tensor-to-scalar ratio interpolate 
from natural inflation to exponential inflation such as R^- and Higgs inflation and brane inflation, 
where the spectral index asymptotes to ns = 1 — 2/N ~ 0.967 for the e-folding number N = 60. We 
also show that a model with the Dedekind eta function gives a sizable running of the spectral index 
due to modulations in the inflaton potential. Such elliptic inflation can be thought of as a specific 
realization of multi-natural inflation, where the inflaton potential consists of multiple sinusoidal 
functions. We also discuss examples in string theory where Jacobi theta functions and the Dedekind 
eta function appear in the inflaton potential. 

I. INTRODUCTION 

The cosmic microwave background (CMB) is isotropic 
to roughly one part in 10®, and its temperature 
anisotropies have the correlations beyond the horizon 
size at the last scattering. This strongly implies that 
our Universe experienced an accelerated cosmic expan¬ 
sion, namely, inflation, at an early stage of the evolu¬ 
tion [1-5]. In the slow-roll inflation paradigm [6, 7], 
the quantum fluctuations of the inflaton become scalar¬ 
mode metric perturbations (i.e. density perturbations), 
while the quantum fluctuations of the graviton become 
tensor-mode metric perturbations [8], inducing primor¬ 
dial i?-mode polarization of the CMB. If the primordial 
B-mode polarization is detected, it will provide a defini¬ 
tive proof of inflation, and determine the inflation scale 
to be around the GUT scale. 

The inflaton field excursion during the last 50 or 60 
e-foldings must be greater than or comparable to the 
Planck scale, Mp ~ 2.4 x 10^® GeV, in order to generate 
large tensor mode perturbations within the reach of the 
current and future CMB experiments. There have been 
proposed many such large-held inflation models. Among 
them, there are two simple and therefore important infla¬ 
tion models: quadratic chaotic inflation [9] and natural 
inflation [10]. The natural inflation consists of a sin¬ 
gle cosine function and it asymptotes to the quadratic 
chaotic inflation in the limit of a large decay constant 
F. The Planck constraint on the tensor-to-scalar ratio r 
reads [11] 

r < 0.11 (95% C.L.), (1) 

and the quadratic chaotic inflation model is on the verge 
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of being excluded, while some of the parameter space for 
the natural inflation {F < 5Mp) has been excluded.^ 

There is one interesting class of inflation models where 
the predicted values of rig and r sits near the best-fit val¬ 
ues of the Planck observation: inflation [3] or Higgs 

inflation [14] (see also Refs. [15-22] for related works). 
Of course this does not preclude any other models lying 
somewhere between the natural inflation and B? infla¬ 
tion. In fact, there are various extensions of the quadratic 
chaotic and natural inflation models, where a variety of 
values of rig and r is realized. For instance, if one con¬ 
siders a polynomial chaotic inflation with the potential 
V{(j)) = -I- K(j)^ -\- X(j)'^ • • •, the value of r can be 

reduced to give a better fit to the Planck data [23-26]. 
Similarly, if there are multiple sinusoidal functions with 
different potential heights and decay constants, a vari¬ 
ety of values of rig and r can be realized [27-29], and 
such extension of the natural inflation is called multi¬ 
natural inflation. It was shown in Ref. [27] that, even 
with two sinusoidal functions, the lower bound on the 
decay constant disappears because a small-field inflation 
(axion hilltop inflation) is realized for a certain choice of 
the parameters. Note that, although multiple sinusoidal 
functions are generated in the extra-natural inflation [30] , 
it does not significantly improve the bound on the decay 
constant as the higher order terms are suppressed. It 
was recently pointed out that this problem can be cir¬ 
cumvented by adding appropriate matter fields in the 
bulk [31], and this provides another UV completion of 
the multi-natural inflation. 

In this paper we propose another extension of the nat¬ 
ural inflation. The inflaton potential in natural inflation 


^ The BICEP2 result [12] prefers a large r which is consistent with 
the quadratic chaotic inflation. However, the subsequent Planck 
polarization data [13] showed that the dust contamination in the 
BICEP2 region is larger than previously thought. 
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is given by a single sinusoidal function, which is a spe¬ 
cific type of periodic function. A periodic potential is 
attractive from the model building point of view because 
it often appears when the flatness of the inflaton poten¬ 
tial is protected by an approximate shift symmetry: the 
inflaton potential which arises from the breaking of a 
continues shift symmetry respects the residual discrete 
one, leading to such periodic potential. We therefore fo¬ 
cus on a class of inflation models with periodic potential. 
Specifically, we consider Jacobi elliptic functions, Jacobi 
theta functions, and the Dedekind eta function as the in¬ 
flaton potential, and we call such inflation models based 
on elliptic functions and other related functions as elliptic 
inflation. 

As an example of the elliptic inflation, we will first con¬ 
sider a Jacobi elliptic function as the inflaton potential: 

F(0) = 2A4(^l-cn2(^A,fc^^, (2) 

where cn(x, k) is the Jacobi elliptic cosine function with 
elliptic modulus 0 < k < 1, and A and F are model 
parameters. We call the inflation model with the above 
potential cnoidal inflation. Interestingly, the cnoidal in¬ 
flation is reduced to the natural inflation in the limit of 
A: —>■ 0, while it is reduced to the exponential inflation 
including the brane inflation, inflation, and the Higgs 
inflation, in the limit of fc —^ 1 . 

The Jacobi elliptic functions are closely related to the 
Jacobi theta functions, which often appear in theories 
with periodic boundary conditions. Periodicity or mod¬ 
ular invariance play important roles in field theories in 
association with symmetries and dualities [32, 33]. In 
string theories, periodicity appears through compactifi- 
cations of extra dimensions when a compact manifold has 
periodic properties. For instance, a holomorphic gauge 
coupling at one loop level is given by the Jacobi theta 
function 'di{v,q) in toroidal compactification [34]. Here, 
n is a open string modulus (diagonal modes of the ad¬ 
joint representation on D-branes), and q = is given 
by a complex structure of a torus r in type HB mod¬ 
els (see also Ref. [35] for their interpretation and fur¬ 
ther generalization). In our later example, the inflaton is 
identified with the real part of v. Another example is a 
Yukawa coupling among matter-like fields, which can be 
described by the Jacobi theta functions 1 ^ 3 ( 1 ;, g) in mag¬ 
netized (or intersecting) D-brane models on tori [36, 37].^ 
The theta functions appear in the gauge and Yukawa 
couplings because of the periodicity of the tori. Simi¬ 
lar Yukawa couplings are obtained in Heterotic orbifolds 
through world-sheet instantons [39].^ Therefore, a scalar 


^ In both IIB and IIA models on tori, it is possible to realize 
magnetized and intersecting D-brane system with any branes if 
one does not require the preservation of supersymmetry [38] . 

® In Ref. [40], they considered a case in which Yukawa couplings 
are given by the Dedekind eta function of Kahler moduli with 
the target space modular invariance. 


potential which depends on the theta functions appears 
naturally, and we shall show that it takes the following 
form, 



( 3 ) 


where a and b denote rational numbers. We call the infla¬ 
tion model with the above potential form theta inflation. 
We Hnd that the theta inflation exhibits similar behavior 
to the cnoidal inflation. 

It is possible to consider inflation models based on 
other special functions. In particular, inflation models 
with the Dedekind eta function was studied in Ref. [41]. 
(See also Ref. [42] for inflation models with modular func¬ 
tions.) Interestingly, the inflaton potential can be ap¬ 
proximated by natural inflation with modulations, lead¬ 
ing to a sizable running of the spectral index. In fact, 
it was shown in Ref. [43] that a sizable running of the 
spectral index can be generated if there are small modu¬ 
lations on the inflaton potential. (See also Refs. [44-50] 
for related work.) Here we will study the inflaton poten¬ 
tial with the eta function in a greater detail. 

The rest of this paper is organized as follows. In Sec. H 
we study phenomenological aspects of the elliptic infla¬ 
tion, where we calculate the predicted values of Ug and r 
for the above three elliptic inflation models. In Sec. HI 
we discuss possible origins of the Jacobi theta functions 
and the Dedekind eta function in the string theories and 
implications for the large-held inflation. The last section 
is devoted to discussion and conclusions. 


II. ELLIPTIC INFLATION 

In this section we Hrst study two examples of the ellip¬ 
tic inflation. One is based on the Jacobi elliptic inflation 
(cnoidal inflation), and the other is based on the Jacobi 
theta function (theta inhation). As we shall see below, 
they provide a one-parameter extension of the natural 
inflation, interpolating from natural inhation (approxi¬ 
mately) to the i?^-inhation. Finally we study the inha¬ 
tion based on the Dedekind eta function, and show that 
a sizable running of the spectral index can be generated 
because of the small modulations on the inhaton poten¬ 
tial. 


A. Cnoidal inflation 

Let us consider the cnoidal inflation model with the 
potential Eq. (2). The Jacobi elliptic functions sn(a;, A) 
and cn(a;, k) can be thought of as an extension of trigono¬ 
metric functions, sin(a:) and cos(a;), respectively. In the 
limit of fc —>■ 0 , they become 

sn(a;,/c) —>• sinx, (4) 

cn{x,k) ^ cosx, (5) 
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FIG. 1. The inflaton potential V(<?!>) with the Jacobi elliptic 
function given by Eq. (2) for k = 0, 0.6, 0.9 and 0.999. 


while, in the limit of fc —^ 1 , 

sn(a:, fc) —>■ tanhx, ( 6 ) 

cn(a:, fc) —>■ secha:. (7) 

In these limits, the inflaton potential V{(j)) becomes 



ns 


{ A^^l —cos-^^ for k ^ 0, 

2A4 tanh^ A ^ 2A4 (l - 

for —>• 1 , 

where we have assumed ^ ^ f in the second equality. 
Therefore, the elliptic inflaton potential Eq. (2) interpo¬ 
lates from natural inflation to the exponential inflation. 
See Fig. 1 for the inflaton potential V {(p) for several val¬ 
ues of k. One can see from the figure that the potential 
becomes flatter around the potential maximum as k ap¬ 
proaches unity. 

To see how the predicted values of the spectral index 
Us and the tensor-to-scalar ratio r change as a function of 
the model parameters, we numerically solve the inflaton 
equation of motion, 

P + SHp + V'iP) =0, (9) 

with V{(j)) given by Eq. (2). The spectral index and the 
tensor-to-scalar ratio can be evaluated by the standard 
formulae. 


Hs ^ 1 + 2r] - 6 e, ( 10 ) 

r ~ 16e, ( 11 ) 


where e and 77 are the slow-roll parameters, 
1 

2{V{P)J ’ 

^ !/(</>) 


( 12 ) 

(13) 


Here and in what follows we have adopted the Planck 
units, Mp = 1, and the prime denotes the derivative 


FIG. 2. The spectral index and the tensor-to-scalar ratio 
in the elliptic inflation Eq. (2) for k = 0, 0.4,0.8,0.99 and 
0.999 from top to bottom with the e-folding number N — 
60. The prediction of i?^-inflation is shown by a red disc. 
We also overlay the Planck result (Planck+WP+HAO) [11] 
in the bottom panel for comparison. The elliptic inflation 
interpolates from natural inflation to i?^-inflation as k varies 
from 0 to 1. 


with respect to the inflaton field (p. In the actual nu¬ 
merical calculation, we have used a refined version of the 
above formulae including up to the second-order slow-roll 
parameters. 

In Fig. 2 we show the calculated values of Us and r for 
several values of k by varying F. One can see that the 
(ns,r) predicted in the elliptic inflation indeed interpo¬ 
lates from natural inflation to exponential inflation. The 
latter includes the R^- or Higgs inflation, in which the po¬ 
tential is given by V{p) ^ (1 — 2 e“V?‘^). Here we have 
fixed the e-folding number to be A^ = 60. We also show 
in Fig. 3 how the Us and r evolve as the decay constant F 
changes. For k sufficiently close to unity (fc > 0.99), the 
spectral index Ug stays around 1 — 2/N = 0.967 for a wide 
range of the decay constant including even sub-Planckian 
values. 


B. Theta inflation 

Next we consider another example of the elliptic in¬ 
flation, the theta inflation with the potential (3). It is 
known that the Jacobi elliptic functions are closely re- 
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FIG. 3. The dependence of ris and r on the decay constant 
F in the cnoidal inflation Eq. (2). 


lated to the Jacobi theta functions, t)- For 

instance, there are following relations: 


where 


k = 


sn(M, k) = 
cn(M, k) = 

^ 2 ( 0 , g) 

^?3(0,g) 


or equivalently. 


V = q = e'-^^, T = 


^3(0,g)t^l(t^,g) 
1?2(0,g) -do('i^,g)’ 
^ 0 ( 0 , g) "92{v,q) 

^?2(0,g) '&o{v,q)" 


u = 7rz;(6<3(0,g))^, 


2K{k) ’ 


K{k) 


(14) 

(15) 

(16) 

(17) 


and K{k) is the elliptic integral of the first kind. The 
Jacobi theta functions are given by 

°° 2 t ,1 

Q / \ \ '' Trim + 27 rm(-u+w ) 

'&o{v,q) = e ^ ’ 

n— — co 

CO 

= l+2^(—l)"g" cos2n7r?;, (18) 

n—1 

n— — co 

CO 

= 2^(-l)"g("+^)%in(2n + l)^t;, (19) 

n—0 


/ \ 2 / \ 

^2(t;,g)= 

n— — co 

CO 

= 2 ^g("+ 2 ) cos (2n + 1)TTu, (20) 

n—0 


d^{v,q)= ^ gVrirn=+2™« 

n— — co 

CO 

= l + 2^g" cos2n7ru. (21) 

n—1 

Here, |g| < 1, i.e., Ini(T) > 0. The Fourier series con¬ 
verge very fast for |g| -C 1 as the higher order terms are 
exponentially suppressed. 

The Jacobi theta functions appear in gauge and 
Yukawa couplings in string theories with toroidal com- 
pactifications. As we shall see in the next section, the 
Jacobi theta functions appear naturally in the scalar po¬ 
tential. 

To be specific, let us consider the following inflaton 
potential, 

m = (22) 

which corresponds to Eq. (3) with a = 0 and 6 = 1. 
Note that a similar inflaton dynamics is obtained for 
the potential with the other theta functions because 
of the identities between the theta functions: '^ 3 ( 1 ' -I- 
1 / 2 , g) = do{v,q), ■d 3 {v + T/2,q) = g), 

'& 2 {v + 1/2, g) = —di{v,q) and so on. In the limit of 
g —> 0 , the above potential is reduced to the natural in¬ 
flation since the higher order terms in the Fourier series 
are suppressed. 

H((()) —>■ ^1 — cosfor g — 7 > 0. (23) 

On the other hand, for |g| —> 1, many sinusoidal functions 
add up so as to make the inflaton potential extremely 
flat around the potential maximum. This is essentially 
same as the axion hilltop inflation where the inflaton po¬ 
tential becomes flat due to cancellation among multiple 
sinusoidal functions [27-29]. In this sense, the theta in¬ 
flation can be thought of as a specific realization of the 
multi-natural inflation. (Also note that the Fourier se¬ 
ries converge fast unless jgj is extremely close to unity.) 
We show the inflaton potential for g = 0, 0.7,0.8 and 0.9 
in Fig. 4. As one can see from the figure, the inflaton 
potential around the local maximum becomes flatter as 
g approaches unity. 

We similarly solve the inflaton dynamics with inflaton 
potential (22), and estimated and r. The results are 
shown in Figs. 5 and 6 . The spectral index Us stays 
around 0.967 for sub-Planckian values of A, if g > 0.7. It 
is interesting to note that the predicted values of Us and 
r become closer to the i?^-inflation as g becomes larger, 
but they do not exactly matches with the i?^-inflation: 
the spectral index tends to be slightly smaller for the 
same value of r. 
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FIG. 4. The inflaton potential V{4>) with the Jacobi theta 
function given by Eq. (22) for 5 = 0, 0.7, 0.8 and 0.9 from 
bottom to top. 



FIG. 5. (ws, r) for the theta inflation with the potential (22). 


C. Inflation with Dedekind eta fnnction 

Now we study an inflation model with the Dedekind 
eta function, 




FIG. 6. Same as Fig. 3 but for the theta inflation (22). 


potential (25) is shown in Fig. 7 for several values of C. 
One can see from the figure that the inflaton potential 
has a maximum at ^ = nF and it receives modulations 
foi' C" < 1. Here, = (0,1) implies the self-dual 

point under the modular transformation of r —> —\/t, 
i.e., C —>■ l/C. For larger values of C, the above poten¬ 
tial asymptotes to the natural inflation with the decay 
constant F ^ as one can see from the expansion of the 
inflaton potential. Note that the normalization of the 
potential height is chosen for visualization purpose. This 
does not affect the predicted values of r, Us and its run¬ 
ning in the following analysis. 


vir)=q^Y[{l-q^n, (24) 

n—1 


where q = and Im(T) > 0. It is known that the 
eta function is given by the special values of the theta 
functions, for instance, ? 7 ^(t) = 5 d 2 ( 0 , r)'d 3 ( 0 , r)do(0, r). 
To be specific, we will study the inflaton potential: 


Vi^) = -A" 


-2 


ttF 


(j) + iC \ + c.c. 


const., (25) 


oc cos - 


F 


+ e cos 4-’ 


The predicted values of Us and r are shown in Fig. 8 
for C = 1.5,1 and 0.85, where we have varied the value 
of the decay constant F, and set the e-folding number 
N to be 60. For (7 = 1 and C = 0.85, the predicted 
values of {ns,r) rotate counterclockwise as F decreases. 
In contrast to the previous two cases, this model does not 
asymptote to the exponential inflation. The dependence 
of Us and r on the decay constant F is shown in Fig. 9, 
and we find that the super-Planckian decay constant is 
required to be consistent with the Planck data. We will 
return to the issue of super-Planckian decay constant at 
the end of this section. 


where (7 is a real parameter, and the constant term is 
added so as to make the inflaton potential vanish at the 
origin, P(0) = 0. In the second line, we have expanded 
the eta function in terms of |g| = e~'^^ <C 1. The inflaton 


Finally let us evaluate the running of the spectral in¬ 
dex, dus/dhik. Recently the inflation model with the 
Dedekind eta function was shown to generate a sizable 
running of the spectral index [41]. The running can be 
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FIG. 7. The inflaton potential with the Dedekind eta func¬ 
tion (25) for C = 1.3,1.1, 0.9 and 0.65 from top to bottom. 



ns 


FIG. 8. (ua, r) for the inflation potential with the Dedekind 
eta function (25). 


expressed in terms of the slow-roll parameters as 


where 


dris 

dink 


-24^2 -b 16er; - 2^, 




(26) 


(27) 


The current constraint on the running by Planck and 
WMAP polarization reads [51] 


dus 
din fc 


=-0.013 ±0.009 (68%). 


(28) 


For e and r] of 0(0.01), the first two terms in Eq. (26) 
are of 0(10“^), one order of magnitude smaller than the 
center value of the current bound. In order to generate 
a sizable running of order 0(0.01), therefore, ^ must be 
large. For a simple class of inflaton potentials expressed 
by Taylor series with the finite truncation, however, a 
large negative running that is more or less constant over 
the observed cosmological scales would quickly terminate 




FIG. 9. Same as Fig. 3 but for the inflaton potential (25). 



F/Mp 


FIG. 10. The running of the spectral index dua/dlnk as a 
function of the decay constant F in the inflation model (25). 


inflation within the e-folding number N <30 [52]. There 
is a loophole in this analysis, and it was pointed out in 
Ref. [43] that a sizable and almost constant running over 
the CMB scales AIVcmb 8 can be generated in an 
inflaton potential with small modulations, because the 
modulations can give a dominant contribution to V'", 
while V, V and V" are not significantly modified. The 
running of the spectral index as a function of the decay 
constant is shown in Fig. 10. For slightly smaller values 
of C, \dns/d\Yik\ can be larger than 0.01, and it is more 
or less constant over the CMB scales. 
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III. POSSIBLE UV COMPLETION 

We have introduced the Jacobi elliptic functions and 
theta functions in the scalar potential to extend the nat¬ 
ural inflation, motivated by the recent Planck data. We 
have shown that they predict the spectral index and the 
tensor-to-scalar ratio interpolating from natural inflation 
to exponential inflation including the B?- or Higgs infla¬ 
tion. Such functions indeed appear in the string the¬ 
ories with extra dimension containing tori. As noted 
already, gauge and Yukawa interactions can be written 
in terms of the Jacobi theta functions in toroidal com- 
pactifications in intersecting or magnetized D-brane sys¬ 
tem [34, 36, 37]. Also on the magnetized orbifolds, e.g., 
T®/(Z 2 X Z 2 ) and x T^/Z 2 , Yukawa couplings are sim¬ 
ilarly expressed with the theta functions [53]. This type 
of function comes from the wave functions of the fields on 
the torus, which respects the periodicity, in magnetized 
brane models (or summation of world-sheet instantons 
in intersecting brane models).^ So, the scalar potential 
with the Jacobi theta function will appear naturally in 
the low energy scales. 

In passing, we note that the theta functions can also 
appear in scattering amplitudes at one-loop level of string 
world-sheet together with the Dedekind eta function (in¬ 
volving moduli integrals) [54]. This is similar in the field 
theories [55]. Also the eta function can appear in one- 
loop threshold corrections to gauge couplings after inte¬ 
grating out heavy Kaluza-Klein modes in the extra di¬ 
mension [56]. Note that the eta function is given by the 
special value of the theta functions, as we have seen be¬ 
fore. 

In the following, we discuss several possibilities to gen¬ 
erate the scalar potential with the theta functions. The 
inflaton candidate is an open string modulus on the D- 
branes, and hence it originates from the Wilson lines, 
positions of D-branes or their linear combination. For 
instance, if a four cycle 5 in a Calabi-Yau, wrapped by 
D7-branes, has the Hodge numbers with 7 ^ 0 or 

/i^’°(iS) 7 ^ 0, there appear Wilson line moduli 

or /i^’°(5) position moduli, respectively [57]. To discuss 
whether such models are really viable for inflation, how¬ 
ever, it is necessary to construct an explicit model of the 
moduli stabilization, supersymmetry (SUSY) breaking, 
as well as the embedding of the standard model sector 
etc. The open string inflation and the related issues have 
been studied extensively in the literature. The purpose of 
this section is to present several possible scenarios where 
the Jacobi theta functions appear in the scalar potential, 
motivated by the recent Planck results on Ug and r. 


The wave functions are invariant under the periodic shift on the 
torus up to gauge transformation, which is represented by double 
pseudo periodicity of the Jacobi theta functions. 


A. Field theoretic cases 

As a toy model, we first consider a case where the 
Yukawa coupling is given by the theta function 'dsiv^q)^ 
where q = , and r is a torus complex structure mod¬ 

ulus r in magnetized brane model ® (or the volume mod¬ 
ulus and B-field in the intersecting brane system), and v 
is identified with the inflaton (j)/2'KF. In a more realistic 
scenario, one has to worry about the stabilization of the 
imaginary component of w, which we shall discuss in the 
next subsection.® 

Now the question is how to generate a scalar potential 
which depends on the Yukawa coupling. One possibility 
is to use a strong dynamics: 

£ = d3{v,q)(p'ilj‘''ip 

y = ^3{v,q){ip){xl}‘'xl}) = (29) 

where (p and ip + are scalar and vector-like fermion 
fields in the strongly coupled hidden sector respectively. 
Here, we have assumed that the scalar develops a non¬ 
zero vacuum expectation value (VEV) and the fermion 
pair condensates through the strong dynamics. 

In a supersymmetric case, we may consider the lYIT 
model [58, 59], where there are four superfields of 
(bJ = 1 ) 2 ,3,4) with fundamental representation 
under SU{2) gauge theory and six singlet superfields 
(j)d = xhe low energy superpotential becomes 

W = (30) 

with Pf(['I'i'I' 7 ]) = Agy^^ 2 )- Here, [^'J'l'j] is a composite 
operator and Asu{ 2 ) is dynamically generated scale under 
the SU{2) theory. Thus, scalar potential is given by the 
F-component of $: 

y — ^su{2) |d3(a,g)|2+0(4>). (31) 

Here, we have assumed that $ = 0 is the minimum ow¬ 
ing to the Kahler potential K = 


^ A Yukawa coupling consists of a product of theta functions if 
a compact extra dimension is given by a product of tori. We 
will focus on the contribution from a single torus and neglect 
contributions from the other tori, assuming that other moduli 
become heavy owing to flux compactiflcation. One may consider 
T'^/Z2 X T2, in which Wilson line moduli could be required to 
explain the smallness of the observable Higgs mass, although the 
moduli are not necessary to explain Yukawa textures [53] . On the 
other hand, is often used for realizing a large (top) Yukawa 
coupling because it contains unity and 1^2 is used for a smaller 
top Yukawas; see Eq. (21). 

® Note also that theta inflation with |g| 1 would require a small 

complex structure, i.e., Im(r) <C 1. Then we may have a diffi¬ 
culty to control the coupling expansion, because then a complex 
structure on a torus can be translated into the volume under T- 
duality. Such a problem may be avoided, if q in the theta function 
depends on a linear combination of moduli or a large volume ex¬ 
tra dimension is realized with stabilized moduli, compensating 
small values of complex structure moduli. 





in which c = 0(10“^) is obtained after integrating out 
heavy modes with the masses of order Agjj(^ 2 )- Alterna¬ 
tively, one may consider the trilinear scalar term of the 
Yukawa interaction, assuming the scalars develop non¬ 
zero VEVs by some dynamics. In this case the scalar po¬ 
tential is approximately given by Y ^ 'm 3 / 2 'd 3 (?;, q) 


B. Supergravity models with moduli stabilization 


In a more realistic stringy set-up one has to consider 
both real and imaginary components of z;. In order to re¬ 
alize an effective single-field inflation, one needs to ensure 
that the other degrees of freedom including the imaginary 
component of v are stabilized with a mass heavier than 
(or comparable to) the Hubble parameter during infla¬ 
tion. To see this, let us consider a supergravity model in 
type IIB orientifold with flux compactification [60-67] . 

The following is the simplified effective action 

K = K{S + S-T + f)-3 log[-z(r - r)] 

+ + ’ (32) 

Oi 

IT = y GAi}{S, r, Va) + A{S, r, 

+IThidden Yukawa (^a ; ^■> '!') + •••, (33) 


where .S' is a dilaton-axion, T is a Kahler modulus, r is 
a complex structure modulus, and Va is a open string 
modulus. Fa is a decay constant for open string mod¬ 
ulus, e.g., = 3/(T -I- T) and Fp = 1/(5' -I- 5) for 

D7 Wilson line moduli and D7 position moduli, respec¬ 
tively [57], and they are expected to be comparable to 
the string/Planck scale. We shall return to the issue of 
the size of the decay constant later. G is the three form 
flux, H is holomorphic three form and the second term 
in the superpotential arises from instantons or gaugino 
condensations. We have included trilinear superpoten¬ 
tial Whidden Yukawa between matter-like fields 'hi in the 
hidden sector. In the Kahler potntial, the ellipsis include 
the terms for and SUSY breaking sector and the quan¬ 
tum corrections [68] in addtion to the Standard Model 
(SM) fields. The ellipsis in the superpotential contains 
contributions from SUSY-breaking and the SM sector. 

With the above effective action, r and 5 can be stabi¬ 
lized by a choice of three form fluxes G: f GAfl D t‘^+St. 
We may obtain a light open string moduli v, a certain 
combination of Va, if 


_9 

dv 


G A fl{S, T, Va) = 0. 


(34) 


The other degrees of freedom in {va} are assumed to be 
stabilized with a heavy mass by / G A D and there¬ 
fore can be integrated out in the low energy. For instance, 
the Wilson line moduli are not included in flux-induced 
superpotential because of the internal gauge symmetry 
[67]. Alternatively, only a part of the position moduli 
on D-branes may be fixed by a selection of the closed 


string fluxes, magnetic (worldvolume) fluxes and geom¬ 
etry through the tadpole condition [69]. (Similar argu¬ 
ment can be applied also to the complex structure mod¬ 
uli [41, 70] and to type IIA moduli stabilization with¬ 
out metric fluxes [71].) Finally, through instantons or 
gaugino condensations on the stack of branes, the Kahler 
moduli T can be stabilized taking account of the SUSY 
breaking. The SUSY breaking also stabilizes the remain¬ 
ing moduli and produces the scalar potential of the Ja¬ 
cobi theta function through a gaugino mass or a trilinear 
scalar term. Here, we expect that Im(z;) obtains the grav- 
itino mass m 3/2 oc f G A fi via the SUSY-breaking effect 
because the Kahler potential depends only on the imag¬ 
inary part of V as in Eq. (33) at the tree level. Thus, 
Re(z;) = (//27rF will be the lightest degrees of freedom, a 
good candidate for the inflaton. 

In the following we consider various ways to generate 
the inflaton potential, which depends on the Jacobi theta 
functions, di{v,q). We will discuss also an issue of one- 
loop breakdown of the shift symmetry along Re(z;). 


Case 1: Gauge coupling depending on the theta function 

Let us hrst focus on the dependence on v through a 
gauge coupling in the strongly coupled hidden sector. 
This n-dependence appears in the non-perturbative term 
as we shall see below. 

On the stack of N D7-branes, the holomorphic gauge 
coupling at one-loop level contains the Jacobi theta func¬ 
tion in the presence of e.g. a stack of n D3-branes [34]: 

Tl 

/one—loop ^ 7^ Q)); 

ZTT 

where q = is given by a torus complex structure 
modulus r and we have assumed that the gauge coupling 
depends on u, and the tree level coupling is given by 
/tree = T. This is regarded as a backreaction from n D3- 
branes to the worldvolume of N D7-branes [35]. In this 
case V denotes the position of D3-branes on the torus. 
When the N D7-branes undergo gaugino condensation^, 


^ Similar type of the superpotential has been observed through 
D-instanton contributions in Ref. [72, 73] with 5L(2,Z) dual¬ 
ity in type IIB supergravity compactified on Calabi-Yau orien¬ 
tifold with odd parity Kahler moduli Q\ For instance, W ~ 
^m{Q(iS) may be generated by D-instantons. Here, 
0m{Q^iS) = theta function of 

weight s/2 with index m. L is some positive definite rational 
lattice of dimension s, and r is some vector which admits an 
expansion in a basis of L with rational coefficients. Then, we 
have dg f G A Q = 0, hence Re(^) can be an inflaton candi¬ 
date if exists, because there will be just Im(^) dependence in 
the Kahler potential at a perturbative level. This looks similar 
to the superpotential discussed in inflation model with Dedekind 
eta function. 
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the second non-perturbative term in Eq. (33) becomes 


^^non—pert OC exp 




(36) 


extra dimension with magnetic fluxes in the string theory. 
The superpotential is given by 


lE = t?3(^;,9)«'r«'j«'fc, (39) 


— ^SU{N) ( <?) 


(37) 


Here, Asu{n) is the dynamically generated scale under 
the pure SU{N) super Yang-Mills theory. Now we have 
the gaugino mass term TO 3 / 2 AA in the SU{N) sector, the 
scalar potential is given by 


y ~ 7713/2( i?i(u,g)) +^{^su{N))^ (38) 


where we neglected the presence of matter d'i and higher 
order terms in Asu{ n) ■ The inflaton is identified with the 
real component of u, the position of D3-branes: Re(u) = 
(j)/2'KF. Thus, the inflaton potential depending on 
arises from the interactions between n D3 branes and N 
D7 branes through gaugino condensation in this case. 

If the inflation is driven by the above inflaton poten¬ 
tial, the Hubble scale during the inflation will be given 
by ^ m3/2Agjjf^j^y In the KKLT case [60], we ex- 
pect < 1713/2 < Asu{n), and so, iJjnf < 1713/2- 

Therefore we may be able to avoid the decompactifica- 
tion problem [74]. The potential with the double gaug¬ 
ino condensations can ameliorate this problem. (See also 
discussions in Ref. [75].) The scalar potential induced 
by a similar mechanism was used to drive slow-roll infla¬ 
tion in the context of the warped brane inflation [35, 76- 
79], where inflation is driven by moving D3-branes in the 
warped geometry with anti D3-branes at the tip. On 
the other hand, here we consider the inflation with a 
D-brane modulus in the bulk extra dimension, and the 
inflaton potential arises from the interactions between 
D3-branes and D7-branes. We may use anti D3-branes 
at the tip of a warped throat to uplift the AdS vacuum 
to the dS one with a small and positive cosmological con¬ 
stant. The interaction between branes and anti branes 
can be neglected if it is a sufficiently strongly warped 
throat. Alternatively, we may consider a field-theoretic 
spontaneous SUSY breaking for uplifting [80-83]. In ei¬ 
ther case the uplifting sector is required to preserve an 
approximate shift symmetry along inflaton direction [76] 
for successful slow roll inflation. 

Lastly, we note that there are u-dependent one-loop 
corrections in the Kahler potential [ 68 ], which induce 
the mass term for the inflaton V ^ 01713/24’'^ with an 

expectation of c = 0(10“^) or so. For successful slow- 
roll inflation, such corrections must be sufficiently small, 

c ^ 10 ^(A|j/(^)/m3/2). 


Case 2: Yukawa coupling depending on the theta function 

As noted already, Yukawa couplings are given by the 
Jacobi theta function in the toroidal compactification of 


where 'I'i,j,fc are hidden matter-like superfields localized 
by the magnetic fluxes on the extra dimension, u is a open 
string modulus and g is a torus complex structure mod¬ 
uli. All the moduli other than v are assumed to be sta¬ 
bilized with a heavy mass. Taking account of the SUSY 
breaking effects, there will be a SUSY-breaking trilinear 
scalar term: UsuSY-breaking = m 3 / 2 i? 3 (?^, <?)^'i^'j^'fc+c.c.. 
(Here, note that are scalar component of the re¬ 

spective superfields.) The scalars can develop non-zero 
VEVs, for instance, via D-term potentials involving with 
Fayet-Iliopoulos term ^ generated by the magnetic fluxes 
and stabilized T: D oc £_ + ftl'l'iPi where qi is an 
U{1) charge of 4'^. See, e.g.. Refs. [64, 84-88] for moduli 
stabilization with D-term potential. Then the potential 
becomes: 

U ~ ad3(u, (?)-|-6 ^i?3(u,q)^ 4-0(4'^). (40) 

Here, Re(v) = (j)/2-KF is a candidate for the inflaton, 
a = 0(7713/24'^), b = 0 ( 47 ^), and we have added SUSY 
F-term contribution Vsusy oc 19<pIUl^, and neglected su¬ 
pergravity contributions of 0(4'®), which contain deriva¬ 
tives of the theta function. Whether decompactification 
occurs depends the VEVs of {4'^}, and a,h < 771^/2 is re- 
quird. As discussed in the previous example, we need 
to suppress one-loop corrections to the inflaton mass 
for successful slow-roll inflation: c < lQ~^a/171^/2 or 
c < 10“^6/m3^2- Here c is a coefficient of the one-loop 
induced inflaton mass term normalized by the gravitino 
mass squared. 


Case 3: Yukawa coupling depending on the theta function 
and strong dynamics 


Now we consider a scenario similar to the case 2, but 
4^^ are stabilized by a strong dynamics as in Refs. [64, 89]: 


Hnon-pert = {Nc — Nf) 


Ae 


-2ttT 


^det(4''=4') 

fkhidden Yukawa — ''^3(77, (?) 4 ?47 4 ^, 


,( 41 ) 

(42) 


where A includes Tli-dependent term in the non- 
perturbative effect as discussed in the case 1, 4''^ -I- 4^ 
are vector-like superfields charged under SU{Nc) gauge 
group with the number of Nf{< Nc) flavors, and 4> is a 
singlet superfields. Note that there are no open string 
moduli (adjoint modes) in the strong coupling brane to 
realize an asymptotically freedom. Here, let us consider 
three stacks of D-branes labeled by (a, ,5, 7 ) for simplic¬ 
ity; 


V = laVa + IpVp + IjVj 


(43) 
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is a linear combination of open string moduli {va, vp,v^) 
on the respective branes [36, 37], and Re(u) = cj)/2TrF is 
the inflaton. Here, la = Ind(-D)/ 3 ^, Ip = Ind(I?)^Q,, and 
= Ind(Zl)a,g, where Ind(Il)a /3 € Z is the index of Dirac 
operator, which depends on magnetic flux between the 
a- and /3-th D-brane and so on®. Hence la = —Ip = Nf 
and I^ counts the number of singlet $. Note that Vj = 0 
when there exists the strongly coupled gauge theory on 
the 7 -th D-brane, because we have assumed there is no 
open string moduli on y-th D-brane. Note that q is given 
hy q = exp(z7rT|/Q,/^/.y|), where r is a complex structure 
moduli on the torus. 

Here, we shall consider stabilization of matter-like 
fields. $ can be fixed by D-term potential in the presence 
of Fayet-Iliopoulos term as discussed in the case 2:® 

D oc e-|d>|2 ~ 0. (44) 

This implies that $ multiplet is eaten by the D(l) gauge 
multiplet on which we are focusing. Here, we have used 
the canonical wave function for simplicity and ^ is com¬ 
parable to the cutoff scale. Thus (<i>) plays a role of 
heavy mass against vector-like pairs, and they condense 

1 27r 

as {'dsiv^q))^ We inte¬ 

grate out them and obtains low energy effective potential 

^ f N 

Wes = N,A^e-^^{Mv,q))^{^)^ (45) 

= ^SU{N,)^^{'^3{v,q))^. (46) 

Here, A 5 ( 7 ( 7 v„) is defined as the dynamically generated 

scale in the pure SU{Nc) Yang-Mills theory in the low 
energy. Thus, low energy potential can be rewritten as 

V ~ q)) through the gaugino mass 

and the trilinear scalar term between $ and 4''^ -|- 
if A does not depend on v. We expect, however, that 
two stacks of flavor branes can give an effect on the 
worldvolume of strong coupling brane, hence we will find 
A ~ [i?i(z;c(,9)]"“[i3i(u/3,g)]”'’, whereg = Ha and 
denote the number of a-th and /3-th D-branes for flavors: 
Two 'di terms come from backreaction of a-th and ,5-th 
D-branes to the y-th brane. Thus the scalar potential is 

V ^ 'm3/2^^SU(Nc) 

Here we have assumed that Va + vp is stabilized by 
closed string fluxes around the origin, v/Nf = Va — vp, 
'di{—v,q) = —'di{v,q), and we took na = np = 1. This 
potential is similar to the case 1, but now it is accompa¬ 
nied with the i^s-dependent term. Note that the decay 
constant is enhanced by 2Nf in the i3i term compared to 
that in the Ds term. . 


f V \ 

Nc 




'&3iv,q) 


® Dirac operator is given by D^p = d + iAa — iAp, where Aa 
denotes the background of gauge potential on the a-th D-brane 
and so on. 

® D will develop the VEV as D ~ 


C. Potential with a gauge coupling depending on 
the Dedekind eta function 


Now let us consider inflation models involving the 
Dedekind eta function, in which a complex structure 
modulus r is an inflaton candidate [41] and a run¬ 
ning spectral index is generated naturally. See also 
Refs. [90, 91] and [92] in Heterotic cases. We choose closed 
sting fluxes such that dr f G A fi = 0 in Eq. (33) with 
assumption that all open string moduli are stabilized. It 
is known that the a holomorphic gauge coupling includes 
the Dedekind eta function [34, 56]: 

/one-loop D log(r7(r)). (48) 

Here this is obtained after integrating out open string 
modes stretched among D-branes which are propagating 
on a part of the bulk extra dimension and hence b denotes 
1-loop coefficient of beta function for the gauge coupling 
in A/” = 2 SUSY sector. We expect that, as in the case 
1, the gauge coupling can depend on di((u),g), where 
q = when there is another stack of n D-branes. 

Including these contributions to the gauge coupling with 
/tree = Y, the gaugino condensation in the pure SU{N) 
theory is given as 


Wnon-pert OC exp “ ^ ^ ^Og{r]{T)) 

= ^SUiN) 




(49) 


Here, A 5 { 7 ( 7 v) is the dynamically generated scale in the 
SU{N). Though the gaugino mass term depending on 
the gravitino mass, the scalar potential includes 


V ~ m 3 / 2 A|y(^) [r]{T)] ^ (z3i((u), q)) ^ -b c.c. 

+ 0(A|^(^)). (50) 


In Refs. [90, 91] on Heterotic otbifolds, b/N = —2 can be 
obtained with a requirement of the modular invariance 
along T direction. If r is a kind of overall modulus, we 
find b/N = —6. 


D. Super-Planckian decay constant 

The decay constant for the inflaton is expected to be 
of order the string/Planck scale. It is, however, possi¬ 
ble to realize a decay constant larger than the Planck 
scale by the alignment mechanism in the presence of 
multiple light moduli [93]. For instance, when there 
are two Yukawa couplings and two light brane moduli 
Vi and V2, i.e., J G A — 0 and Whjdden Yukawa — 

i?3(7iUi -I-/ 2 U 2 , q)'i>^ Filsilivi + T 2 U 2 , g)4''®, an enhanced 
decay constant can be obtained with JJi — J 2 I// < 1, 
Ii = I2 = I and (T') < (T). The inflaton becomes 
Re(z;i — V2) and the decay constant is enhanced by a fac¬ 
tor of // j Ji — J 2 1. Indeed, we have seen the enhancement 
of the decay constant by 2Nf in the scalar potential in 
the above case 3. 
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IV. DISCUSSION AND CONCLUSIONS 

The Universe must be reheated by the inflaton decay 
some time after inflation. The inflaton appears in the 
gauge and Yukawa couplings, and so, it may decay into 
the SM gauge bosons (or gauginios if kinematically acces¬ 
sible) through the gauge couplings or into the SM matter 
and the Higgs or right-handed neutrinos and B — L Higgs 
through Yukawa couplings. For instance, the following 
operators will be relevant to the reheating process: 

^y^^B-LNN. (51) 

Here, we have taken {(j)) ^ 0 because the modulus devel¬ 
ops a non-zero VEV in general, and the decay constant 
Uu is comparable to the string scale, which is not nec¬ 
essarily equal to the decay constant in the inflaton po¬ 
tential. t, H, ifiB-L and N denote the SM gauge 

bosons, top quark, B — L Higgs and right-handed neu¬ 
trino respectively. The thermal leptogenesis is viable if 
the reheating temperature is enough high as 10® GeV 
[94]. The last operator will be useful to implement non- 
thermal leptogenesis, even if the reheating temperature 
becomes lower [95, 96]. 

In this paper we have studied an extension of the nat¬ 
ural inflation, where the inflaton potential has a peri¬ 
odic property. We have proposed elliptic inflation where 
the inflaton potential is based on the elliptic functions 
and related functions. Focusing on the elliptic and theta 
inflation, we have shown that the predicted Hs and r 
interpolate from natural inflation to the exponential in¬ 
flation including i?®-inflation, where Us asymptotes to 
ris ~ 1 — 2/N ~ 0.967 for N = 60. We have also studied 
the inflation model with the Dedekind eta function and 


shown that a sizable running of the spectral index can 
be generated owing to the small modulations on the in¬ 
flaton potential. For a certain value of the parameter, it 
is possible to realize dus/dlnk ^ —0.01 which is more or 
less constant over the CMB scales. 

Periodicity often appears in the field theories and the 
string theories through compactifications. Interestingly, 
such periodic backgrounds restrict the allowed form of 
the potential, and as a result, the inflaton potential takes 
a specific form rather than simple cosine function(s). For 
instance, toroidal background forces infinite series of co¬ 
sine functions to add up to give the Jacobi theta functions 
(cf. Eqs. (18) - (21)). Indeed, the Jacobi theta functions 
are known to appear in gauge and Yukawa couplings in 
toroidal compactifications, where the variables are open 
string and complex structure moduli. We have explored 
various ways to generate the scalar potential which de¬ 
pends on the Jacobi theta functions. The inflaton is iden¬ 
tified with the real component of the lightest open string 
modulus, and we have discussed the stabilization of the 
other moduli fields based on a supergravity model in type 
HB orientifold with flux compactification. 
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